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Abstract. The presence of strongly correlated superconductivity in M117S15 has 
been recently established. In this work, we compare the normal and superconducting 
parameters of a single crystal of RI117S15 with those of a polycrystalline sample of 
PdirSeis which is reported here for the first time to be a superconductor below 2.2 K. 
PdivSeis, which is iso-structural to Rh^Sis (space group Pm3m), has very different 
properties and provides for an interesting study in contrast with Rh^Sis. We see that 
large unit volume of Pd^Seis and the large separation of Pd-Pd atoms in its structure 
as compared those in a pure Pd metal lead to the absence of strong correlations in this 
compound. Finally, we report band structure calculations on these two compounds and 
compare the density of states with estimates from heat capacity data. Upper critical 
field, Heat capacity (down to 500 mK), Hall effect and band structure studies suggest 
that RI117S15 is a multiband superconductor. 



PACS numbers: 74.70.Ad, 74.25.Bt, 74.25.Dw, 74.25.F, 72.80.Ga, 71.27.+a, 71.28. +d, 
71.20.Be, 71.15.Mb, 71.15.Ap 
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1. Introduction 

Recently there has been intense activity in the study of transition metal sulphides 
and selenides due to a variety of their properties such as, strongly correlated 
superconductivity, spin and charge density waves. In our effort to understand the 
physical properties of noble metal chalcogenides, we have studied cubic RJ117S15 
(Ref. [U |2]) which is considered as an excellent catalyst by chemists (Ref. [3j 0]) and 
Pdi 7 Sei 5 . 

Rh 17 Si5 (known as Miassite), a mineral mainly found near the Miass river in Russia, 
and PdiySeis (known as Palladseite), a mineral mainly found in Brazil, are iso-structural 
cubic compounds whose crystal structure is shown in fig. [TJ They have a large unit cell 
with 64 atoms in each (two formula units) and Rh (or Pd) occur in four inequivalent 
crystallographic sites. Some of us had reported strongly correlated superconductivity 
(T c = 5.4 K) in Rh 17 Si5 (Ref. p]) with magnetization, resistivity, heat capacity and 
Hall measurements. The compound was remarkable for its enhanced susceptibility, 
large Sommerfeld coefficient (108 mJ/mol-K) and unusually high upper critical field 
(H C 2) values (around 13 T at 3 K) well above the Pauli limiting field (9.9 T). We had 
ascribed the origin of strong correlations to the possible occurrence of a narrow Rh 
d-band of states at the Fermi level possibly due to the close lying Rh atoms at the 3d 
and 6e positions. The distance between these atoms (2.58 A) is even lesser than that 
between the closest Rh atoms in Rh metal (2.69 A). In iso-structural Pdi7Sei 5 , however, 
the lowest Pd-Pd distance (2.78 A) is more than that in elemental Pd (2.75 A). Pdi7Sei5 
is found to be a superconductor below 2.2 K. 

Our paper is arranged as follows. The first section deals with the sample growth 
and structural characterisation details. The second has results and discussion on the 
superconducting properties of Rh 17 Si 5 . The third deals with the superconducting 
properties of Pd^Seis which are compared with those of RI117S15. In the fourth section 
we report band structure calculations on these systems. The fifth section contains the 
conclusions. 

2. The Samples 

The compounds were prepared by heating stoichiometric compositions of the two 
elements at the rate of 8 deg/hr to a temperature of 1150 deg C in an alumina crucible 
that was sealed in an evacuated quartz tube. Thus reacted mixture is annealed at 1050 
deg C for two days and then cooled down firstly at 8 deg/hr till 600 deg C and further 
at 30 deg/hr. The samples were characterised using powder XRD and the single crystal 
of RI117S15 showed clear cubic spots of Pm3m space group in a Laue diffractometer. The 
lattice constants were estimated to be around 9.9093(2) A for Rh 17 Si 5 after a Reitveld 
analysis. 

Since this is the first report on superconductivity in Pdi^Seis, we have confirmed 
that the compound is indeed in a single phase by both powder XRD and Electron Probe 
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Micro Analysis (EPMA). In fig. [2] we show the refined XRD pattern for Pd^Se^. 
There is no impurity peak and the pattern refines to the cubic Pm3m structure of 
PdiySeis with a lattice constant of 10.6044(7) A by a Reitveld analysis. Hence, the 
XRD pattern indicates an essentially single phase sample. Phase sensitive EPMA shows 
almost uniform stoichiometry with Pd:Se concentration of 52.25:47.75 atomic percent 
in the analysed region (an area of 1.5mm x 2mm). A 17:15 ratio of Pd and Se should 
give a 53.1:46.9 atomic percent ratio. Providing for a small error in estimation, we 
conclude that the majority of the sample has indeed formed in the desired 17:15 ratio. 
However, small regions of Pd:Se concentration of 49.62:50.38 atomic percent are seen. 
Inspite of small off-stoichiometric regions, our results indicate that there is indeed bulk 
superconductivity in Pdi 7 Sei 5 . 

The specific heat capacity measurements were performed in a Physical property 
measurement system (Quantum Design, USA). The magnetisation measurements 
were done in a Magnetic property measurement system (Quantum Design, USA). 
Transport measurements were performed in a home-made system and in an Adiabatic 
demagnetisation refrigerator (Cambridge magnetic refrigeraton, UK). 

3. Superconducting Properties of single crystalline Rh^Sis 

3.1. Specific Heat Capacity measurements 

In fig. [3] we show the zero field specific heat capacity (C p ) data with C p /T plotted 
against T 2 . The continuous line is a fit to 

C p = 7 T + (3T 3 + 5T 5 (1) 

in the low temperature region in the normal state where 7 is the electronic 
contribution, is the lattice contribution and S is the contribution due to anharmonicity. 
We obtain a 7 of 108.41 mJ/mol-K 2 , (3 of 0.657 mJ/mol-K 4 and 5 of 0.00139 mJ/mol- 
K 6 . The presence of a significant anharmonicity could be a result of a complex 
phonon structure in these compounds which could arise largely due to a complex crystal 
structure. Using an equal entropy estimation, we can estimate a T c of 5.36 K. From j3 
we obtain a Debye temperature {9d) of 456 K which indicates a hard lattice. Using Or, 
and McMillan's equation (Ref. [5]) with a Coulomb parameter (fi* = 0.13) we obtain an 
electron-phonon coupling constant (A) of 0.58 which puts Rh 17 Si5 i n the intermediate 
coupling regime. However, from the jump at the transition (around 1558.554 mJ/mol- 
K), ACp/7 T c is around 2.68 which indicates that the electron-electron interaction 
is in the strong-coupling regime as opposed to a weak-coupling BCS value of 1.43. 
Therefore, we conclude that phonon-mediated mechanism is insufficient to account for 
strong electron correlations in Rh 17 Si5. Using the value of 7 one can also compute the 
enhanced density of states N*(0) from N*(0) = 0.4248 7 states/eV-formula unit where 7 
is in mJ/mol-K 2 . With the earlier stated 7 value we get N*(0) = 46 states/eV-formula 
unit. 
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The inset in fig. [3] shows very low temperature heat capacity data with C p /T 
plotted against T. Clearly the low temperature values of C p /T seem to level off at a 
residual non-zero value (7 res ~ 16 mJ/mol-K 2 ) whereas we expect the heat capacity to 
go to zero as the temperature approaches OK. This behaviour has been seen in multi- 
band superconductors like MgB 2 and more recently in the pnictides (Ref. [19]). Such 
systems typically have multiple gaps in the superconducting state where some gap is 
small enough to contribute to the heat capacity even at low temperatures. 

3.2. Upper Critical Field measurements 

Fig. [4] shows the upper critical field of Rh 17 Si 5 as a function of temperature from 70 
mK to 5.4 K as estimated from ac and dc magnetization measurements. The very low 
temperature data (below 2 K) have been included from a recent study (Ref. j6]). From 
this one can estimate the value of the upper critical field at H c2 at K to be 20 T which 
is in agreement with the recent single crystal study of Settai et al (Ref. [7]). The slope 
(dH c2 /dT) near the superconducting transition is estimated to be - 3.6 T/K which is 
smaller than the value quoted by Settai et al (Ref. [7]). From H c2 (0), we can estimate 
the BCS coherence length (£ ) using H c2 (0) = (fi /2irt;l which gives us a £ of 4*10~ 7 
cm. 

A strong-coupling BCS theory estimation of the orbital limit of the upper critical 
field of a type II superconductor is the WHH theory expression (Ref. [8]) 

H c2 (0) = -0.69T C 

from which we get H c2 (0) to be 13.3 T which is much smaller than the actual value. 

The Pauli paramagnetic limiting field (Hp(0)), which refers to the energy 
scale required to depair the Cooper-pairs, is roughly 1.84*T C (in Tesla) for a BCS 
superconductor. This value is around 9.9 T for Rh^Sis which is half the actual value 
of H c2 (0). Taking into account the electron phonon coupling (A) we can estimate an 
enhanced Pauli limit (Ref. [9]) H P (0)* = 1.84T C (1+A) 1/2 which gives a value of 12.4 
T for RI117S15 which is again much smaller than the actual value. For RI117S15 to be 
Pauli limited in a strong coupling BCS sense, the coupling parameter A will have to 
be much larger than 3. It can be speculated that the reason for such a large upper 
critical field could be strong electron-electron correlations (as in the heavy fermions) 
or unconventional pairing or large spin-orbit coupling (as in the Chevrel Phases) or a 
combination of the three. We cannot conclusively decide on any one reason based on 
the available data. We have also estimated the penetration depth Ao as 7000 A from the 
/iSR measurements performed at PSI, Switzerland recently (Ref. [TO])- Using Ao and £0 
we estimate the Ginzburg Landau parameter k to be around 175 which makes Rh 17 Si5 
an extreme type II superconductor. From this one can estimate a lower critical field of 
16 Oe using H cl (0) = folnn/^l. 

We also note that the positive curvature of the H c2 curve near T c is typical of multi- 
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band superconductors like MgB 2 (Ref. [TT]) and the Borocarbides (Ref. [12J). Using the 
emperical expression 

H c2 = H* c2 (l - T/T c ) 1+a 3 

for quantifying this positive curvature (Ref. [HI [12]), we get a good fit from 
0.5T C <T<T C as seen in fig. H] with a H* 2 of 30.3 T and 'a' of 0.24 (a value similar 
to Borocarbides and MgB2 in clean limit). This, the low temperature heat capacity 
data (previous subsection) and the fact that the Hall data in Rh 17 Si 5 shows a sign 
change in majority carriers (Ref. [TJ) strongly hint at a multi band scenario for Rh 17 Sis. 

4. Superconducting properties of polycrystalline Pd 17 Se 15 

4-1. Resistivity measurements 

Fig. |5] shows the temperature dependence of the resistivity (p(T)) of Pdi 7 Sei5 from 1.6 
to 300 K. Inset (a) shows the transition region (T c = 2.2 K for Pdi 7 Sei 5 ). Inset (b) 
shows a T 2 fit to the low temperature data (p = Po + A*T 2 ) which gives a value of 43.2 
fiQcm for the residual resistivity (po) and a value of 0.014 pf2 cm/K 2 for A. The value 
of A for Pd 17 Se 15 is comparable to that of Rh 17 Si 5 Ref. [I]. 

The lower graph in fig. [5] shows the resistivity data of the single crystal of Rhi 7 Si5 
for comparision. The two curves are distinctly different in the appearance of a knee-like 
feature at around 60 K in Rhi 7 Sis which is completely absent in Pdi 7 Seis which has a 
conventional metallic behaviour. We do not, at present, understand the reason for the 
knee-like feature. However we note that this kind of behaviour has been observed in 
many transition metal based superconductors, particularly in many A- 15 compounds 
(Ref. [13]) which have a sharp feature in density of states at Fermi level. In that 
scenario, the Fermi level moves to drastically lower values of density of states at higher 
temperatures and causes a knee in the resistivity. Such a mechanism could be in 
operation in Rh 17 Si5. We also note that the sign change in majority carriers as reported 
in (Ref. pQ) occurs at around 60 K as well. 

4-2. Specific Heat Capacity measurements 

Fig. [6] shows a plot of C p versus T for the polycrystalline Pdi 7 Seis in the absence of 
applied magnetic field from 2 to 15 K. Our range of C p measurement, unfortunately 
cannot capture the transition region well. The solid line in the same figure is a fit to 
equation (1) described earlier. The fit yields a 7 of 22 mJ/mol-K 2 , (3 of 4.9 mJ/mol- 
K 4 and of 0.11 pJ/mol-K 6 . From the value of /3 we obtain a Debye temperature (9d) 
of 233.3 K which indicates that Pd 17 Sei5 is a much softer lattice as compared with 
Rh 17 Si5. From the values of T c and T c , we estimate a value of 0.54 for the electron- 
phonon coupling constant (Ref. [5]) for Pdi 7 Seis which is very similar to the value of 
Rh 17 Si5. The anharmonicity is smaller than that of Rh 17 Si5. More importantly, the 
enhanced density of states, N*(0), as calculated from the earlier expression, is around 
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9.35 states/eV-formula unit and is smaller than that of Rh 17 Si5 by a factor of 5 indicating 
that the correlations are much weaker in Pd 17 Se 15 . The reduction in the density of 
states in PdiySeis as compared to that of RI117S15 is also reflected in the temperature 
dependence of normal state magnetic susceptibility as shown in the fig. [7J 

4-3. DC magnetisation measurements 

Fig. [7J shows the dc susceptibility curves of RI117S15 and Pd^Se^ at different 
fields. Since none of Rh, Pd, S or Se carry a magnetic moment we expect to see 
a temperature independent susceptibility, essentially sum of Pauli paramagnetism, 
Landau diamagnetism and core diamagnetism. However, although the values of 
susceptibility (x(T)) are small, the data display a distinct temperature dependence for 
Rh 17 Si 5 . There is a gradual increase in susceptibility as one goes to lower temperatures 
until the superconducting transition. The weak temperature dependence of x(T) could 
arise due to a sharp density of states at the Fermi level which leads to the temperature 
dependence of the Pauli spin susceptibility. We believe that the mechanism is similar 
to the observed in the A-15 compound, VsSi (Ref. [14J). This conjecture is recently 
supported by the observation of temperature dependent Knight shift from the 103Rh- 
NMR experiment Ref. [15]. However the x(T) of PdiySeis is diamagnetic with impurities 
(presumably at few ppm levels) dominating at low fields. The value of x(T) is nearly two 
orders smaller than that of Rh 17 Si 5 indicating that the density states is much smaller 
which is confirmed by the heat capacity data. 

4-4- Upper Critical Field measurements 

Finally, we show the temperature dependence of the upper critical field (H^T)) of 
PdiySei5 from 270 mK to 2.2 K in Fig. |HJ The curvature of H C 2 nearT c is more like a 
conventional superconductor as opposed the one seen in Rh 17 Si5. However, both show 
a tendency to saturate at lower temperatures. The saturation is better seen in RI117S15 
because we have data upto 70 mK whereas, we have data only upto 270 mK in PdiySeis. 
The H C 2(0) is estimated to be 3.3 T by an extrapolation which is nearly 6 times smaller 
than that of RI117S15. The BCS coherence length £ is 100 A as estimated from the 
H c2 (0). The slope (dH c2 /dT) near T c is -1.5 T/K which yields a value of 2.3 T for 
H C 2(0) from WHH theory (equation (2)). However, the Pauli paramagnetic limiting 
field (Hp(0)) is estimated to have a value of 4 T which suggests that that orbital critical 
field is Pauli limited in contrast to the situation in RI117S15. The upper critical field 
data indicate the Pd^Se^ is a conventional superconductor as opposed to Rh 17 Si 5 . 

5. Band Structure Calculations 

In this section, we present some band structure calculations on these two compounds. 
Firstly, we will present a LDA (Local Density Approximation) calculation of the density 
of states and band structure and compare with the experimental data. Later we will 
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present an LDA+U (which takes into account Coulombic interaction between electrons) 
calculation and compare that with experiments. 

5.1. LDA calculations 

These calculations were performed using the WIEN2K code (Ref. [16]) which is a Full 
Potential Linearised Augmented Plane Wave calculation (FPLAPW). Calculations were 
also performed with a Full Potential Local Orbital code (FPLO) (Ref. [TT1 ITS] ) and 
the results are very similar to the WIEN2K results. Also the above calculations were 
performed in an LDA framework as well as in a Generalised Gradient Approximation 
framework and the results match very closely. We will be presenting the WIEN2K 
calculations in an LDA framework here. 

In fig. [9] we have the band structure of Rh 17 Si5 juxtaposed with the Density of 
States plot. We see that the Fermi level (here E^ = 0) indeed lies on a peak of the 
DOS. We can see both electron-like and hole-like pockets at the Fermi Level supporting 
the fact that there are both types of charge carriers in the system. This also supports 
the conjecture that RI117S15 could have a strong compensation due to oppositely charged 
carriers as speculated by Settai et al from the magnetoresistance data (Ref. [7]). Around 
the M direction, we also see flat electron-like and hole-like bands which could contribute 
to the peak in the DOS at Fermi level. Such flat bands also hint at a multi-band scenario 
supporting the Hall data and the fact that the upper critical field curve has a multi-band 
like positive curvature. 

Fig. HO] shows the DOS plot in the total energy range of calculation. The DOS 
at Fermi level is around 27 states/eV- formula unit. We recall that the ehanced DOS 
from the heat capacity data was around 46 states/eV- formula unit. From this one can 
estimate the electron-phonon coupling (A) using 

^ , DOS from heat capacity . ^ ^ 
DOS from band structure 

which gives a A of 0.7. This value is slightly larger than that calculated from 
Mcmillan's formula (0.58). 

In fig. [10] we have also indicated the contributions by Rhodium d-bands and 
Sulphur p-bands to the DOS. Nearly 80 percent of the total DOS is contributed by the 
Rh d bands of which the largest contributors are the Rh24m atoms mainly because of 
the large number of such atoms. Interestingly, when one calculates a contribution per 
atom, it turns out that the Rh24m atoms still contribute more than twice the other 
Rh or S atoms. This is in contrast with our earlier speculation that it must be the 
closer lying Rh3d and Rh6e atoms that contribute to the DOS at Fermi level. However, 
this supports the NMR data which says that that the Rh 24m atoms contribute the 
maximum towards the Knight shift (Ref. [T5]). 

In fig. CD] we have the band structure of Pdi7Sei5 juxtaposed with the DOS plot. 
We see that the Fermi level doesn't lie on a peak unlike Rh 17 Si5. There is also, unlike 
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Rh 17 Si5, an absence of flat bands at Fermi level which which could be why the Fermi 
level doesn't lie on a peak of DOS. 

Fig. [12] has the DOS plot in the entire range of calcualation. Here, as in RI117S15, 
the Pd d bands contribute more than 70 percent of the DOS at Fermi level. However, 
though the Pd 24m atoms are the major contributors to the total DOS, the Pd lb atom's 
contribution per atom is larger than all the other Pd and Se atoms. We see that the DOS 
at Fermi level is around 16 states/eV- formula unit. This is a remarkably large value 
in comparison with that from the heat capacity data (which is around 9.3 states/eV- 
formula unit). This is a puzzling discrepancy since the enhanced DOS from heat capacity 
data is always expected to be larger than that obtained from an LDA calculation. To 
check if strong correlations are the origin of this discrepancy we performed LDA + U 
calculations on both the systems. 

5.2. LDA + U calculations 

When strong correlations were taken into account in the calculations, the DOS at Fermi 
level decrease for both the compounds. The Fermi level still lies on a peak for RI117S15 
and in a valley for Pdi7Sei 5 . But the value of DOS at Fermi level is around 6 states/eV- 
formula unit for Pdi 7 Seis and 11.4 states/eV-formula unit for Rh^Sis. This, when 
compared with the heat capacity estimates gives a A of 0.56 for Pdi 7 Seis and 4 for 
Rh 17 Si5. While the A seems to match a McMillan's equation like estimate for Pdi 7 Sei5, 
it seems unrealistically large for Rhi 7 Sis. We comment on these calculations in the next 
section. 

6. Conclusions 

Earlier (Ref. [1]), we had conjectured that the strong correlations, as deduced from the 
Sommerfeld coefficient, enhanced susceptibility and large upper critical field, in Rhi 7 Si5 
were due to strong Rh-Rh interactions which form a large density of states at the Fermi 
level. On the other hand Pdi 7 Sei5, whose superconductivity we report for the first time 
here, has much weaker electron correlations as deduced from smaller values of the same 
parameters as mentioned above. Since the unit cell volume of Pdi 7 Seis is larger than 
Rhi 7 Si5, the interactions between Pd atoms could be weaker here and hence the weaker 
correlations. 

The deductions from the band structure calculations, however, are not 
straightforward. It appears that an LDA calculation seems to suffice to estimate the bare 
density of states of Rhi 7 Sis while an LDA+U calculation seems necessary for Pdi 7 Seis to 
match the A values got from McMillian's equation. On the other hand LDA for Pdi 7 Sei5 
and LDA+U for Rh 17 Sis seems to compare very badly with A from McMillian's equation. 
We believe that McMillian's equation, which is a based on a strong coupling BCS theory, 
is a reasonable estimate for Pd^Seis since it is a conventional superconductor. Whereas 
McMillian's equation may not be applicable to Rhi 7 Sis since it is an unconventional 
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superconductor. This was also seen when the coupling parameter derived from the 
jump at the transition was much larger than A from McMillian's equation. As we 
saw from the upper critical field, it would require a A much larger than 3 to be Pauli 
limited in the BCS sense which is unrealistic. We believe that an LDA+U calculation 
is necessary to estimate the DOS of these systems and this gives a good estimate for 
Pdi 7 Sei 5 . However, DOS from an LDA+U calculation for Rh 17 Si 5 and its comparision 
with the DOS from heat capacity, cannot give a A matching a Mcmillan's equation 
estimate since Rh^Sis is a strongly correlated system. It will interesting to study other 
families of chalcogenide superconductors with narrow 4d and 5d bands and compare for 
similarities with the '17-15' family, (Rhi 7 Si 5 and Pdi 7 Sei 5 ). 
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Figure 1. [colour online] The crystal structure of RI117S15 and PdiySeis contains two 
formula units per unit cell 
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Figure 2. [colour online] The powder XRD pattern of PdiySeis as refined with the 
calculated Pm3m structure of Pdi7Sei5 by a Reitveld analysis 
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Figure 3. [colour online] C p /T versus T 2 for RI117S15. The continuous line is a fit to 
the equation mentioned in the text. Inset shows the low temperature C p /T versus T 
data which levels off at 7 res ~ 16 mJ/mol-K 2 




Figure 4. [colour online] H C 2 versus temperature for RhiySis. The inset shows 
the positive curvature at the transition. The continuous line is a fit to the equation 
mentioned in the text 
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Figure 5. [colour online] Resistivity data of Pdi7Seis and RI117S15. Inset (a) shows 
the transition region and inset (b) shows a T 2 fit. The lower figure is the RI117S15 
curve showing the knee-like feature 
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Figure 6. [colour online] C p /T versus T 2 for Pdi7Sei5. The continuous line is a fit to 
the equation mentioned in the text 
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Figure 7. [colour online] DC susceptibility curves of Rh^Sis an d Pd^Seis at different 
fields. Note the two different scales for the two compounds 
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Figure 8. [colour online] H C 2 versus temperature for Pd^Se 
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Figure 9. [colour online] The band structure of Rh^Sis. The total DOS is indicated 
on the right side of the plot 
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Figure 10. [colour online] DOS for RIT17S15 from LDA calculations. Contributions of 
Rh and S are indicated in red and blue respectively 
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Figure 11. [colour online] The band structure of PdiySeis. The total DOS is indicated 
on the right side of the plot 




Figure 12. [colour online] DOS for Pdi7Seis from LDA calculations. Contributions 
of Pd and Se are indicated in red and blue respectively 



